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Abstract 

Kaluza-Klein photon in universal extra dimension models is one of the most attractive dark 
matter candidates as a weakly interacting massive particle. Having a characteristic split spectrum 
in split universal extra dimension the relic density of Kaluza-Klein photon with 900GeV mass is 
in good agreement with the observed dark matter amount in our Universe. Interestingly Kaluza- 
Klein photon in the same mass range also provides excellent fits to the recently observed excesses 
in cosmic electron and positron fluxes. The amount of gamma-ray contributions, mostly from tau 
decays, can be significant around 300 GeV, thus can be observed or constrained by the forthcoming 


Fermi-LAT diffuse gamma-ray data. 
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I. INTRODUCTION 


Dark matter problem is one of the most pressing problems in particle physics and astro- 
physics. Even though the presence of dark matter is firmly established by several indepen- 
dent experiments and observations, we still do not know what is the major component of 
dark matter. 

An attractive dark matter candidate is weakly interacting massive particle (WIMP) as 
it naturally reconciles with the standard thermal history of big bang cosmology. For being 
a proper WIMP candidate a particle is required to be a stable one (or its lifetime has to 
be much longer than the age of Universe) otherwise it could not contribute to the structure 
formation of Universe and its decay products should have been detected. Often a discrete 
symmetry is introduced to make a WIMP stable: R-parity in supersymmetric models, T- 
parity in Little Higgs models and Kaluza-Klein (KK) parity in extra dimension models. Once 
a discrete symmetry, e.g. Zə, is conserved, the lightest Z odd particle is automatically stable 
and can not decay to the standard model (SM) particles which are set to be even under the 
Zə transformation. 

Recently a particular interest has been attracted to Kaluza-Klein dark matter in univeral 
extra dimension model (UED) hl pi ld pi J ld pi ld ld hid Ail and its variety split-UED 
adl 15| since it has been shown that a TeV Kaluza-Klein dark matter might be responsible 
for the excesses in positron fraction, flux of electron plus positron and gamma-ray in cosmic- 
ray data lids from PAMELA Lod Lash, ATIC Ri PPB/BETS ba, HESS and 
Fermi-LAT . Indeed a TeV scale Kaluza-Klein photon, B4, is naturally expected as a 
dark matter candidtate in UED and it has advantages over a Majorana spinor dark matter 
(e.g. neutralino in MSSM) since no helicity suppression arises in annihilation into two light 
fermions. 

Split-UED has been suggested to address the problem of excessive hadron production 
in conventional UED models by invoking a 5D bulk Dirac mass term for quarks which is 
generically allowed by symmetries of the theory in a way of keeping Kaluza-Klein parity 
intact Production of hadron is suppressed by heavy Kaluza-Klein quarks so that the 
better fit to the PAMELA antiproton data is provided in split-UED. 

In this paper we further consider Kaluza-Klein dark matter model with particular em- 


phasis on new electron (here, electron means electron and positron since Fermi-LAT doesn’t 
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tell them apart.) and diffuse gamma-ray data of Fermi-LAT. Three issues are particularly 
interesting here. First, the smoother and further extending shape of e~e* excess in high 
energy regime seems to suggest a rather heavier dark matter (> 900 GeV) than what is 
suggested by ATIC/PPB-BETS (~ 620 GeV). Second, lack of peaky tip of flux also sug- 
gests the direct production of e~e* channel sub-leading. Third, a sizable high energy diffuse 
gamma-ray flux is expected in our model from inverse Compton scattering of electrons as 
well as tau decays. 

This paper is organized as follows. After a brief review of the model in Sec. [I| we first 
consider the relic density of Kaluza-Klein photon in split-UED with splitting spectra of 
leptons and quarks in Sec. Concerning possible experimental constraints when we split 
fermions, we present bounds on 5D bulk masses for leptons in Sec. Then we fit the 
positron fraction and the excess in e~e* flux by introducing a larger Kaluza-Klein scale 
(Sec. IIA). The diffuse gamma-ray excess in high energy regime (2, 100 GeV) is expected 
in split-UED, which could be seen by forthcoming Fermi-LAT diffuse gamma- ray detection 
(Sec. IIB). Then we conclude. 


II. SPLIT-UED (REVIEW) 


Split-UED is an extra dimension model where a natural candidate of dark matter arises 
as a consequence of the symmetry of the theory, dubbed Kaluza-Klien parity. All the 
standard model fields are universally propagating through the bulk so that their Kaluza- 
Klein excited states provide phenomenological implications if the compactification scale is 
accessible by the experiments (1/R ~ TeV). In the minimal setup, the gauge group is 
exactly the same as the one in the standard model, G = SU(3), x SU(2)w x U(1)y, so KK 
gluons, weak gauge bosons and photons are involved. Quarks and leptons are also alleviated 
to higher dimensional ones and their KK spectra are doubled. One chirality of their zero 
modes are projected out by the orbifold condition, leaving the other one as the standard 
model fermions. The most prominent feature of split-UED is the presence of (double) kink 
masses for fermions while keeping the KK parity. Once 5D bulk mass is introduced, Kaluza- 
Klein modes get additional mass contributions and become heavier while zero mode remain 
massless because of the orbifold condition. Let us be more concrete. The five dimensional 


action of split-UED is the same as the minimal UED (MUED) except the 5D bulk mass 
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terms: 


L 
Ssplit-UED = SmUED — fef vg mg (y) UW; (1) 
=Ñ 


where W; are the bulk Dirac spinors containing quarks and leptons in their zero modes and 


the kink mass, m4 (y) = y/0(y), is introduced with a step function defined as @(y > 0) = 1 


and @(y < 0) = —1. One should notice that the zero mode wave functions, ~ et Jms(y)dy, 
are even functions under the inversion about the middle point of extra dimension, y = 0, so 
that Kaluza-Klein parity is respected. Here + signs are determined by chirality. 


Klauza-Klein states of fermions, on the other hand, get additional masses: 
Ma = Mo + ky + H5 (2) 


where the first term (m2) comes from the ordinary SM Yukawa interaction, the second term 
(k2) from the momentum of the extra dimension and the last term (u2) from the 5D bulk 
mass. kn is determined by ps; = +k, cot k„L for KK modes. Again + sign depends on 
chirality. All the details of Kaluza-Klein decomposition are in the Appendix. 

Having the generic idea of split- UED, we can control KK spectra by turning on some bulk 
mass parameters. In order to suppress the hadronic annihilation cross section and avoid the 
unnecessary flavor problems, we first choose our 5D bulk masses for quarks to be universal 
and larger than the typical KK scale, which is pid = [5 0;;. For the charged leptons, in order 
to control their annihilation cross section ratio among different flavors which are dominantly 
coming from the right-handed components, we choose separate 5D bulk masses Heg, Hup and 
Urg to achieve that. For the the left-handed leptons, we assume that their 5D bulk mass 
are generally small so that their couplings are still almost KK conserving. Therefore their 


contributions to the lepton flavor violation and four fermion operators are negligible. 


A. Relic abundance of Kaluza-Klein DM 


Relic abundance of Kaluza-Klein dark matter has been extensively studied (See l as a 
review on Kaluza-Klein dark matter and references therein). It is pointed out that due to 
coannihilation with the right-handed charged leptons (mainly e®), the lightest KK particle 
(LKP) in 500 GeV ~ 700 GeV mass range is needed to generate the right relic density in the 
original MUED U pi i. In split-UED, however, due to the presence of 5D bulk masses 
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for the leptons (even when it is very small), the degeneracies between the LKP and the KK 
leptons are automatically removed. In this case, a LKP mass between 900 GeV to 1 TeV 
will predict the right relic abundance pi Ha from 5 years WMAP data and is exactly 
the necessary value needed to fit the Fermi and Hess data. 

It is possible to even get a heavier dark matter to account for the relic abundance in a 
slightly extended model. For instance, by turning on the brane kinetic term for the SU (2) 
gauge bosons , it is possible that the SU(2) KK gauge bosons have degenerate mass 
spectra with LKP. In this case, the coannihilation effects from the SU(2) KK gauge bosons 
will increase the required masses to account for the right relic abundance up to 1.5 TeV [6]. 

One may wonder why the coannihilation with KK singlet leptons will decrease the dark 
matter mass for the right relic abundance while the coannihilation with SU(2) KK gauge 
bosons will increase it. The reason is that the overall effective cross section, which roughly 


depends on 1/m?,,,, goes like 


N 
gigj 

Oeff = Dij REE , (3) 
ij Jeff 


when coannihilation is included. When one particle species is removed from the degenerate 
spectra, not only we remove the relevant coannihilation channels which will decreases o¢,, 
but we also decouple some particle species which will decrease geff and increases Ceff, SO 
the overall effect is really a competition. The former wins in the coannihilation with the 
SU (2) KK gauge bosons (one has to decrease the dark matter mass to compensate that so 
their relics density is fixed) while the latter wins in the coannihilation with the KK singlet 


leptons due to their small couplings. 


B. Splitting leptons 


In this section, we consider that leptons also have their own 5D bulk masses in order to 
split their KK spectra and the relevant constrains. Because KK number is also violated in 
the lepton sector, the even KK gauge bosons will induce the four fermion operators between 
leptons and between leptons and quarks at tree level. The contact interactions between 
electrons and the SM fermions are parameterized by an effective Lagrangian, Leff 

2 
Leff = aa >. pee Tr (4) 


Ba i, j=L,R 


5 
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where g? is taken to be 47 as the p meson coupling. 6 = 1(0) for f =e (f #e), AG. is the 
scale of the contact interactions, € = +1, e; and f; are left or right-handed spinors. 

Because all the even KK Z bosons are almost completely in the weak gauge eigenstate 
(even KK W? boson), they only couple to the left-handed fermions, while the even KK 
photon couple to both the left-handed and right-handed fermions. So contact interactions 
involves right-handed fermions are less constrained since only even KK photon will con- 
tribute. Hence we only split the right handed electron KK spectrum to maximally suppress 
the electron annihilation branching ratio. 

In our setup, we only introduce the 5D bulk masses to localize of the SM fields at the 
center of the extra dimension 1. All the couplings between the SM fermons and even KK 
gauge bosons have the same signs therefore the interference is constructive except the case 
with right handed up quark (wz) and the quark doublet (Qz) due to the opposite sign of 
hypercharges of them to hypercharge of the right-handed electron (ep). 

We first consider the constrains for the 5D bulk mass of the right-handed electron. If 
the up or/and the down quark has 5D bulk mass, then their contact interactions with the 
electron will give the most stringent bound. From the global fit of the relevant data bd 
when the SU(2); gauge symmetry is assumed, by summing over the contributions from all 
even KK photons, we obtain the constrained Pd tone space for pe, and 1/R in Fig. 


from most stringent bound A, + = 15.2 TeV [29] from the eg-dpr contact interaction 


Soho ten Rofo (de R)[— 2] < SE (5) 
7 " 7 ý 5 3M pon (P) 7 (Avra)? 


The terms of g$ o and g2? o are the ratios of couplings between KK gauge bosons and SM 


fermions to the SM coupling for er and dp, respectively, 
L i 
J- dy fo fo faem fe, = S-i dy fen fy fae 
L > J2n00 — 
I dy fof. fen 


Here we use mgen(R) ~ 2n/R and go) > V2 by assuming all the KK quarks are 


e — 
J2n00 = 


(6) 


-= l 
Ji dy f yo Ly FB 


decoupled. As we can see in Fig. [I] the bound is quite loose and one can completely split 
the electron if the KK scale 


1 
F > 770 GeV. (7) 


1 Tn our previous publications id, we adopted the case of the SM fields localized toward boundaries in 


which there is a sign difference in our calculation. 
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FIG. 1: In the figure on the left we show the allowed space in (peg R,1/R) parameter space. The 
dark region is excluded by bounds for the contact |ererdrdp] interaction. When 1/R > 770 GeV 
basically any value of Hep is allowed. The figure on the right shows the (Heg, Hug) parameter space 


with the KK scales 1/R = 700,800 GeV and 900 GeV where the dark region is excluded. 


Clearly, this bound is satisfied for the dark matter mass used to fit the electron spectra 
in this paper, mzrgp > 900 GeV. 

The constraints on the 5D bulk mass of right-handed muon and tau are mainly coming 
from LEP [80]. The bound for tau is relatively weak, and allows to introduce its arbitrary 


5D bulk mass. The bound for 5D bulk mass of muon can be obtained from 


2 
X ER g 4T 
Fie Hex a el) r] 5 (8) 


B2n ~ (Art)? ” 
by taking Afp, = 11.9 TeV ?. The allowed parameter space on fup vesus different epp is 
presented in Fig[]when we fix the KK scale 1/R = 700,800 GeV and 900 GeV, respectively. 
From the above discussions, we conclude that by introducing the 5D bulk masses for the 
right handed components of electrons and taus, we can adjust the their Bı annihilation 


branching ratio arbitrarily, which can be used to fit the Fermi-LAT electron and gamma-ray 


data. 


2 Compare to Akr = 7.0 TeV, ARr+ = 8.2 TeV, the bound from muon is the most restrictive one bo. 


201708.00321v1 


chinaXiv 


II. COSMIC-RAY FROM KK DARK MATTER 


The KK dark matter, Bı, mainly annihilates into fermion pairs among which charged 
lepton pair is dominant. Stable particles, like proton (p), electron (e~) and photon (7), 
and their antiparticles are generated from the casecade decay of hadrons and heavy charged 
leptons in the B,B, annihilation process and propagate to the Earth. In this section, we 
calculate the electron, positron and gamma-ray from the annihilation of KK dark matter 
and compare with the recent observations from PAMELA and Fermi-LAT collaborations. 
In the framework of split-UED, antiproton flux can be naturally suppressed and agrees well 
with the PAMELA data l, as seen in Ref. h, therefore, we will not show the result of 


antiproton here. 


A. Electron and Positron 


The calculation of flux of electrons and positrons is exactly the same as our previous 


work |15]. Here we only give a brief review and some important equations in our analysis 


for completeness. The propagation of electrons and positrons (e*) can be discribed by a 


diffusion equation as 


K(B) 9 fa (E, P) + LBE) fee (E, P+ QEP) =0, (9) 


where K(E) is the diffusion coeffieient, fe+(E,T) is the density of e* per unit kinematic 


energy, b(E) is the rate of energy loss and Q(E, T) is the source of e+. For annihilation of 


BB, : 
(EF) = Br (22) Den (AP) (10) 


where dN,+/dE is the energy spectrum of e* obtained by using a Monte Carlo program, 
PYTHIA 


6], and the index 7 runs over all quark and charged lepton pairs, and p(7’) is 
the dark matter profile. In our numerical calculations, we adopt an overall boost factor 
Br. There are known sources of boost factor: local clumps in dark matter profile fal ld, 
Sommerfeld enhancement effect by a long range attractive force and the Breit-Wigner 
type resonance effect [34]. In split-UED, without assuming a new attractive force or further 
tuning of mass spectra for large resonance effect, boost factor may mainly comes from local 


clumps. 
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FIG. 2: The positron fraction (et /(e* + e~)) with mpg, = 900 GeV and M1 propagation model, 
compared with PAMELA data. The mass of KK quark is taken to be about 3mg,. The black 
curve is for e: w: 7 =1:1:1 case and the blue curve is for 0.3 : 1 : 0.3 case. . 

The isothermal halo model is adopted as 

m Po 
alo") = ; 11 
Phalo(T) 1+ (r/r)? (11) 

where r = |r] is the distance from our Galactic center, re = 3.5 kpc and po is the parameter 
that is adjusted to yield a dark matter local halo density of 0.3 GeV/cm? [35] in our solar 
system. 

Finally, the fluxes of electron and positron observed near the Earth is given as 


D (E) = TflE, ro), (12) 


where c is the speed of light and rẹ = 8.5 kpc is the distance from Solar system to the 
Galactic center. 
For the background of electron and positron, we adopt approximations Bi of the 


fluxes of primiary electron, secondary electron and secondary positron: 


a glg E E E 
or) = EEN 4 3.952 GeV “cm “sec” ‘sr, 
a 0.7E%7 1 2 ies 
P(E) = Ty THOE + GO0E? papa COV cm sec sr, 
Ap Eo 
(8) = oe _—: 1G lem sec Isp (13) 


1+ 650E?3 + 1500E4? 
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where F is in units of GeV. We also use a free parameter k multiplying the primiary electron 
backgound and we take k = 0.7 in our numerical study. In Fig. 2] and Fig. [B] we show the 
fracton of positron (et /(e* +e7)) and the total flux of electron and positron (et + e7) from 
split-UED model with mg, = 900 GeV and compare with the experimental data, using so- 
called M1 propagation model Bi The mass of KK quark is split to be three times heavier 
than the Bı dark matter through this paper below. The black solid line is for univeral KK 


leptons, namely all of the KK leptons have the same mass and are only slightly heavier than 


Bı. Therefor there exist equal amount of e*, w* and T= in the final state. 

The case of e : 4:7 = 1:1:1 fits the Fermi-LAT data acceptably well with Xot <1 

0l. Moreover, in the set-up of split-UED model, the spectrum of KK leptons can be also 
split by turning on the bulk mass terms, and the fraction of different flavor of SM leptons 
in the final state will be different. Therefore, we can fit the Fermi LAT data even better 
if the rato ise: u : T = 0.3 : 1: 0.3 (in such case, the mass spectrum of KK lepton is 
7) N el) > 1 © Bı) with the blue solid line in Fig. B| It can be seen that the peaky shape 
of the curve is much smoothen and fit the data very well. With such a parameter set, we 
can also fit the PAMELA data in Fig. [2] In conclusion with a 900 GeV KK dark matter, the 
anomalies of positron fraction and total flux of electron and positron can be well explained 


in the split-UED model. 


B. Gamma-ray 


Now, we consider the gamma-ray from the inner galactic plane (IGP, 0.25° < |b| < 4.75°, 
0.25° < £ < 29.75°, 330.25° < Z < 360°) and the intermediate galactic latitudes (’10-20’, 
10° < |b] < 20°, 0° < £ < 360°), both of which are observed _by Fermi-LAT and the 
preliminary data have been presented in several conferences lad , using the parameter 


sets that explain electron and positron data. The gamma-ray can be produced directly 


through the decay of 7° generated in the decay process of r} and the hadrons in the final 
state of B,B, annihilation. The latter is subdominant and can be neglected when the KK 
quark is much heavier than Bı, which is the case for satisfying the antiproton data of 


PAMELA. The flux per unit energy of these primary gamma-ray is given as 
dJ, L J dN, IRT 

Ee = ——— < ov >; | — e)de, 14 

(F _ An ami, >, = ( dE ) f, # l ) ( ) 


10 
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FIG. 3: The total flux of positron and electron with mp, = 900 GeV and M1 propagation model, 
compared with Fermi-LAT and HESS data. The mass of KK quark is taken to be about 3mp,. 


The black curve is for e : u: T =1:1:1 case and the blue curve is for 0.3: 1 : 0.3 case. . 


=> 


where i denotes the channel which produces photon in the final state and f, p7( )d@ is the 
integration along the line of sight (los). 

Other than the primary gamma-ray mentioned above, the Inverse Compton Scattering 
(ICS) in the regions we are considering can also contribute to the energetic gamma-ray since 
there exists hard electrons and positrons in Bı dark matter annihilation. And indeed, we 
will see later that the ICS can be comparable with the primary one. For ICS calculations, 
we follow the procedures in Ref. and a nice review of ICS can be found in Ref. [44]. The 
differential flux of the scattered photon with a energy F is given as 


dJ; 1 < ov al ak P(E ET) ee 
SERA Sy l\dé dE ——_——_Y (E’), 15 
( a) Fie J röd f on) (15) 


where P(E, E’,1) is the differential power emitted to a photon of energy E by a positron or 
electron with energy E’, €(E’) is the total energy loss rate for a positron or an electron with 
energy E’ and Y(E”) is the total number of positron or electron with energy larger than F’. 


The analytical expressions of P and € are 


11 
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0:25" < |b| < 4.75°,|10.0|1.7 x 10711 7.0 x 1075 


0.25° < 29.75°, 


Ls 
330.25° < £ < 360° 


10° < |b| < 20°, | 2.3 |8.9 x 10713 1.3 x 1075 1 


0° < £ < 360° 


TABLE I: The geometrical factor J and parameters of modeling the Interstellar Radiation Field 


in Eq. (19). 


P(E, E',r) = (16) 
ae i A ( 7 ey, a E pam ae eee 7 J 
AE (17) 
of foe go 


where op = 8rr2/3 = 0.6625 barn is the total Thomson cross section, E = E/E’ with y = 
E!/me, n(e,r) is the number density of background photon with energy €, q = me E /(4ey(1— 
E)) and a = 4ey/m,. Approximately, for the region of Galaxy we are interested in, the 
spatial dependence of Eq. (15) can be integrated separately and the simplified result reads 


oo aye Lol jan f pF) a, (18) 
dE J ics E An mh, fis €(E’) 


where J is the geometrical factor, AQ is the solid angle of the observed region. The distri- 


bution of photon bath can be approximated by using the blackbody-like spectra , 


e 1 
e) = > Naip LT (19) 


where the index a represents the spacial region we are interested in and 27 is for the different 
component of the photon bath, i.e. star light (SL), the infrared radiation (IR) from the 
galactic dusts which absorpt the star light and re-radiate out photons and photons of cosmic 
microwave background (CMB).The numerical values for IGP and ’10-20’ are listed in Table [I] 
and the temparture T; for SL, IR and CMB are 0.3 eV, 3.5 meV and 2.725 K, respectively ; 


12 
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FIG. 4: The gamma-ray signal from Bı dark matter annihilation for the regions of Inner Galactic 


plane (upper) and intermediate galactic latitudes (lower). The black solid and dashed curves are 


the sum of primary and ICS contributions when the final state charged lepton branching ratios are 


(0.3 : 1 : 0.3) and (1 : 1 : 1) respectively. The blue and red lines are the gamma-ray from 7* and 


hadrons from Bı dark matter annihilation, respectively, and the magenta lines are contributions 


of three components of ICS in the (0.3 : 1 : 0.3) case. 


In Fig. [4] we show the gamma-ray signal for Inner Galactic plane region and the region of 


intermediate galactic latitudes in the case that the KK leptons are not universal which has 


a better shape in flux of electron and positron shown in Fig. Bland compare with the Fermi- 


LAT preliminary data. The KK quark mass spectra are taken to be Mma ¥ 3/R so that we 


13 
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are able to reach their KK quark production and decay signals at the LHC hs. The black 


curve is the sum of primary and ICS contributions, the blue and red lines are the gamma-ray 


from T7 and hadrons from B, dark matter annihilation, respectively, and the magenta lines 
are contributions of three components of ICS. In the region of low energy (Æ < 10 GeV), 
the dark matter signal is much smaller than the observed data, however, starting from few 
tens of GeV, the signal is about a factor of 2 ~ 3 smaller. A characteristic of our model 


is that a bump at Æ ~ 300 GeV is predicted, due to the main contributions from primary 


T~ and also the ICS from the star light in a subleading way, which can be checked soon if 
higher energy of gamma ray data is available. We also show the case of universa KK lepton 


in dashed line for reference. The behavior is similar to the non-universal case, but with more 


energetic gamma-ray due to more hard e* and more T7. 


IV. CONCLUSION 


We consider the pair annihilation of the lightest Kaluza-Klein photons in split-UED as 
a primary source of recently observed cosmic ray positron and gamma in PAMELA and 
Fermi-LAT. Leptophilic property of dark matter suggested by the PAMELA antiproton 
data is naturally realized in split-UED. 

As the mass of dark matter particle around 900 GeV and its primary annihilation channel 
being lepton pairs with e : u : 7 = 0.3: 1: 0.3 (or 1:1: 1) we successfully fit the all the 
cosmic ray data. A particularly interesting prediction of our model is that the excess of 
cosmic gamma-ray flux, if observed by the forthcoming data of Fermi-LAT diffuse gamma- 
ray, peaks at Æ ~ 300 GeV range. If there is no excess in the high energy region, then 
Fermi-LAT will put an upper bound on the tau fraction in our model. 

Finally we point out another interesting prediction for the collider phenomenology. In 
the case of splitting right handed charged leptons (i.e. 0.3 : 1 : 0.3 case) a large cross section 
of dilepton (in particular, eger and TrTR) production is expected at the LHC through 2nd 
KK gauge boson exchanges. As these leptonic signals are rather clean we expect that the 


detection would be promising and we leave it for future study. 
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Appendix: KK decomposition in split-UED 


Let us consider a massive fermion on an orbifold $1/Z, with two fixed points at y = —L 
and y = +L: 
+L 7 E 
S= fef dy [iY T” 3Y — ms(y) WV] . (20) 
=; 


A y-dependent kink-mass is introduced for keeping KK-parity as 


ms(y) = u 0(y), (21) 


where 6(—L < y < 0) = —1 and 0(0 < y < L) = +1 and gamma matrices are’ = (4#, iy5). 
Left (Right)-chiral fermion is defined as usual ys z/R = FY z/g and a generic Dirac fermion 
is decomposed by Y = WV; + Wp. 

Varying the action with respect to Yz and Vp we obtain the standard bulk equation of 


motion which are given by 


iy ð YL — YUR = M5ÜR = 0, (22) 
iq O UR — YFL = ms Yg = 0 (23) 


then using yY r/r = +V zr we finally get 
(FÔR = ms) Uri + iy"O.V LR = 0. (24) 


Now we would like to discuss how to perform the Kaluza—Klein decomposition of these 
fields. In general, when the fermion belongs to a complex representation of the symmetry 
group, the KK modes can only acquire Dirac masses and the KK decomposition is of the 


form 
Wir =J nainu) (25) 
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where y? jp are 4D spinors which satisfy Dirac equations: 
iq" ô VEJR = Mn YRL (26) 


Plugging this expansion into the bulk equations we get the following set of coupled first 


order differential equations for the wave functions fp IR: 


(+05 — ms) frj + Mnf tyr = 0. (27) 


Applying (Fô; + ms) on the first order equations we get a decoupled second order equa- 
tions in the bulk: 


0 = (F8; + ms) |(F3s — ms) faye + Mn fiyr] 
= (ð — m5 + m, F ms) fay 


= (05+ At) frit, (28) 


where A? = m? — m? F ms. 


1. Zero mode solution: mo = 0 


For n = 0, we can find a massless solution (mp = 0) rather easily. In the bulk (y 4 0), 


the equations in 27] are reduced to a simple first order equations: 


(+05 — ms) fry, = 0 (29) 


having simple solutions 


Pry ly) ~ oF Pam Way — #7 (y) = Nayre Pe (30) 
where the normalization factors are obtained by the normalization condition f an ee JL Fads 


Eu 
1 — eF2HE 


Nryjt = (31) 


Depending on the sign of u the shape of wave functions are determined. If u > 0, e.g., 


f$ is localized toward the middle point (y = 0) and f? towards the end points (y = +L). 
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2. KK mode solution: Heavy modes (m2 > p?) 


Depending on the sign of A? the wave functions fp/,(y) will be either sines and cosines 
or sinhes and coshes. Here we first consider the case with A? = k? > 0. In this case the KK 
modes are heavier than the bulk mass since m? = u? +k? > p?. We call them heavy modes. 


The wave equation for heavy modes looks simple as: 


(35 + ka) frL = 0, (32) 
and their generic solutions are 
Frc y) = Ary COS kny + Spyz sin kny. (33) 
a’s and ĝ’s are related by Eq. 27 
Fah phn — M587 p+ MnBRip = 9, (34) 
tT RKn = MALJR + MnARJL =) (35) 


From the continuity condition at the middle point (lim. .9 ( fiir) — fry nl+€)) = 0) we 


get a useful formula: 


tkn (btr — Brie) 
QLJR 


where we have used a? = a< = a from the continuity condition. 

Now let us consider boundary conditions. We can have two independent choices of Dirich- 
let boundary conditions according to the Z orbifold condition: fg(L) = 0 = fr(—L) (DL) 
or fr(L) =0 = fr(—L) (DR). 


0 = frr(L) = AER cos kn L + Brrr sin kn L, (37) 
0 = frjr(-L) = ATR cos kn L — Bin sin kn L, (38) 
or 
Nn, > _ Qn < 
A O S (39) 
201 jr 


Now combining the continuity condition in Eq. B6]and the Dirichlet boundary condition in 


Eq. BI] we get the master equation: 
p= Fkn cot(kn L) 
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for (DL/DR), respectively. This equation determines the KK spectrum for heavy modes 


with any given values of ju. 
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